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1. Introduction 

Overwhelming evidence supports the existence of dark matter (DM) based on its gravitational 
influence on ordinary matter. Although the DM paradigm is well-established, the nature of the 
particle that may constitute DM is not yet known. One of the leading candidates for the DM 
particle are weakly interacting massive particles (WIMPs) that existed in thermal equilibrium with 
standard model particles in the early universe. The relic density of WIMPs is set by “freeze-out” 
when the rate of WIMP interactions drops below the expansion rate of the Universe. For particles 
with an annihilation cross section on the weak interaction scale such as WIMPs, this process can 
naturally produce a relic density equal to the density of DM observed today. 

Self-annihilation of WIMPs in regions with high DM density can produce stable secondary 
particles detectable by ground- and space-based observatories. The gamma rays produced in these 
annihilations are an especially attractive target for experimental searches. Unlike charged secon¬ 
daries, gamma rays can be traced back to their point of origin enabling DM signals to be more 
easily discriminated from astrophysical foregrounds. The Large Area Telescope onboard the Fermi 
gamma-ray observatory is sensitive in the energy range from 20 MeV to more than 300 GeV and 
is well-suited to search for the gamma-ray annihilation signature of DM. With its excellent angu¬ 
lar resolution, the Fermi-LAT can search for signals in a variety of targets including the Galactic 
Center, dwarf spheroidal galaxies (dSphs), and galaxy clusters. 

Milky Way (MW) dSphs are especially promising targets due to their proximity and the ab¬ 
sence of intrinsic sources of gamma-ray emission. The MW dSphs are also predominantly found 
in high-latitude regions of the sky where the diffuse gamma-ray foregrounds are lower and more 
easily discriminated from a potential DM signal. The DM content of these objects can be robustly 
measured by kinematic modeling of the stellar velocities of member stars. The expected DM signal 
is proportional to the line-of-sight integral of the DM distribution known as the J-factor. J-factors 
for known dSphs generally fall in the range between lO^^GeV^cm^^ and 10^°GeV^cm^^. By 
comparison, models for the DM distribution of the Galactic Center region predict J-factors at least 
an order of magnitude larger. However searches for DM signals in the Galactic Center are compli¬ 
cated by the systematic uncertainties associated with the diffuse astrophysical foregrounds which 
are brightest in this region of the Galaxy. Because the astrophysical foregrounds are much lower 
and more easily modeled in dSphs, DM searches in dSph galaxies can offer comparable sensitivity 
to those in the Galactic Center. 

In this paper we present an updated search for DM signals using the sample of 25 dSphs from 
[1] based on 6 years of LAT data processed with Pass 8. ^ Following the approach of previous 
works [3, 4, 5, 1, 6, 7], we use a joint analysis to combine the data from multiple dSphs and increase 
the sensitivity of our search. We discuss the implications of these measurements in the context of 
the Galactic center excess and other experimental constraints on the DM annihilation cross section. 

2. Pass 8 

One of the main improvements of the present analysis with respect to prior works is the use 

*The content of this proceeding is a condensed version of [2]. See this paper and its supplementary materials for 
additional details. 
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of the Pass 8 data set. Pass 8 is a comprehensive revision to the LAT event-level analysis that 
reconstructs and classifies individual events from the raw data collected by the LAT subsystems. 
With respect to the current Pass 7 Reprocessed data release, Pass 8 contains a number of 
improvements to the event reconstruction including algorithms to identify and remove instrumental 
pile-up, an improved energy reconstruction that extends the LAT energy range below 100MeV and 
above 1 TeV, and a more robust pattern-recognition algorithm for track reconstruction [8]. 

The event classification in Pass 8 has been fully redesigned and provides a completely new 
set of event classes for high-level science analysis. The P a s s 8 event class selections have been 
optimized using boosted decision trees to maximize the rejection power for cosmic-ray background 
events while maintaining a high efficiency for gamma rays. The Pass 8 event analysis also ex¬ 
pands the existing classification framework by adding event types, subdivisions of the event classes 
based on event-by-event uncertainties in the directional and energy measurements. The PSF event 
types are ordered by the quality of the angular resolution from PSFO (worst) to PSF3 (best). At 
3.16GeV the 68% (95%) containment radii of the incidence-angle-averaged PSF for the best and 
worst PSF event types (PSF3 and PSFO) is 0.17° (0.35°) and 0.92° (2.3°). Our maximum likelihood 
analysis of the dSphs combines the four P8R2_SOURCE_V6 PSF event types in a joint likelihood 
function. We estimate that splitting the event sample by PSF event type improves the sensitivity 
of our analysis by ~ 10%. Considering both the improvements in instrument performance from 
Pass 8 and the expanded six-year data set, we expect an improvement of between 1.7 and 2.2 in 
the median sensitivity to WIMP annihilations as compared to the four-year analysis of [1] based on 
Pass 7 Reprocessed data. 

3. Analysis 

We analyzed LAT data collected during the first six years of the mission (2008-08-04 to 2014- 
08-05). From this data set we select events in the energy range between 500MeV and 500GeV 
in the Pass 8 SOURCE class. We apply a zenith angle cut of < 100° to remove photons from 
the Earth Eimb. Time periods around bright gamma-ray bursts and solar flares are also excised 
following the prescription of the Fermf-EAT third source catalog (3EGE) [9]. Data are binned into 
10° X10° square regions of interest (ROIs) in Galactic coordinates centered at the position of each 
dSph. We use a spatial binning of 0.1° and 24 logarithmic energy bins (8 bins per energy decade). 

We perform a standard binned likelihood of the 25 dSphs in our sample using the Eermi Sci¬ 
ence Tools and the P8R2_SOURCE_V6 IREs. We use the standard templates for isotropic and 
galactic diffuse emission. ^ Point-sources within each ROI are taken from the 3PGE catalog [9]. 
Each ROI is first fit with a baseline model in which the normalizations of the diffuse components 
and catalog sources within the ROI are left free to vary. Eor each dSph we use a spatial profile cor¬ 
responding to an NEW DM density profile with halo scale radii and characteristic densities taken 
from [1]. 

After fitting the baseline model to each ROI, we derive a set of bin-by-bin likelihoods following 
the methodology of [1]. In each energy bin, we scan the normalization of the DM source while 

^Specifically we use the v06 release of the Galactic interstellar emission model and isotropic templates documented 
at http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html. 
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leaving the parameters of the background model fixed. The set of bin-by-bin likelihood functions 
for each target is then used to reconstruct the global likelihood for a variety of DM spectral models. 

Given the LAT likelihood function for a dSph, we construct a likelihood for the parameters 
of the DM model (/I) from the product of LAT and J-factor likelihoods. We define the likelihood 
function for target i as 


= {a/,7/} I ^i) I 17obs,i, Oi). 


(3.1) 


where 6/ is the set of nuisance parameters that includes both parameters from the LAT analysis (a/) 
and the dSph J-factor (7/) and ^/ is the gamma-ray data. The second term in Eq. 3.1 is the J-factor 
likelihood function which we use to model the uncertainty in the J-factor arising from the finite 
precision of the stellar kinematic analysis. We use a log-normal parameterization of the J-factor 
likelihood given by 


-^(7/1 7obs,i) 


_i_ X g^(*°glQ(-^')^^°glo(Abs,i))V2CT} 

In (10)7obs,i \( 27r(j/ 


(3.2) 


where 7/ is the true value of the J-factor and 7obs,i is the measured J-factor with error a/. We obtain 
the J-factor likelihood parameterization by fitting a log-normal function with peak value 7obs,i to 
the posterior distribution for each J-factor as derived by [10]. 

We perform a combined analysis for DM signals in our dSph sample by forming a joint like¬ 
lihood from the individual dSph likelihoods. Our combined sample consists of 15 dSphs that have 
kinematically determined J-factors and avoid ROI overlap and corresponds to the same set that 
was used in [1]. The combined likelihood function is given by the product of the individual dSph 
likelihoods, 

^(jl, {0/1 1 &) = I ^/), (3.3) 

i 

where the index i runs over the 15 dSphs in our combined sample. With the combined dSph 
likelihood, we test the hypothesis of a DM signal in the entire sample where the expected signal 
in each dSph is weighted in proportion to its J-factor. For a given mass and annihilation channel, 
upper limits on the annihilation cross section are evaluated with the delta-log-likelihood technique, 
requiring a change in the profile log-likelihood of 2.71/2 from its maximum for a 95% CL upper 
limit [11]. 


4. Results 

We find no evidence for significant gamma-ray emission from any of the individual dSphs or 
in the combined sample. Using the joint likelihood for the 15 dSphs in the combined sample, we 
derive limits on the WIMP annihilation cross section as a function of WIMP mass. Figure 1 shows 
the combined limits evaluated for the bb and annihilation channels. The expectation bands 
show the expected range for these limits based on an analysis of 300 sets of randomly selected 
blank fields with \b\> 30°. We find that the limits from the six-year Pass 8 analysis improve on 
the limits of [1] by a factor between 3 and 5. This improvement can be attributed to a number 
of factors including the improved sensitivity of the Pass 8 data set and the different statistical 
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realizations of the two data sets. Because the Pass 8 six-year and Pass 7 Reprocessed four- 
year event samples have a shared fraction of only 20^0%, the two analyses are nearly statistically 
independent. For masses below 100 GeV, the upper limits of [1] were near the 95% upper bound 
of the expected sensitivity band while the limits in the present analysis are within one standard 
deviation of the median expectation value. 




Figure 1: Constraints on the DM annihilation cross section at 95% CL for the bb (left) and (right) 
channels derived from a combined analysis of 15 dSphs. Bands for the expected sensitivity are calculated by 
repeating the same analysis on 300 randomly selected sets of high-Galactic-latitude blank fields in the LAT 
data. The dashed line shows the median expected sensitivity while the bands represent the 68% and 95% 
quantiles. For each set of random locations, nominal J-factors are randomized in accord with their measure¬ 
ment uncertainties. The solid blue curve shows the limits derived from a previous analysis of four years of 
Pass 7 Reprocessed data and the same sample of 15 dSphs [1]. The dashed gray curve corresponds to 
the thermal relic cross section from [12]. 




Figure 2: Constraints on the DM annihilation cross at 95% CL section for the bb (left) and (right) 
channels derived from the combined analysis of 15 dSphs with 6 years of Pass 8 data. For comparison limits 
from previously published searches are shown from LAT analysis of the Milky Way halo (3(7 limit) [13], 112 
hours of observations of the Galactic Center with H.E.S.S. [14], and 157.9 hours of observations of Segue 1 
with MAGIC [15]. Pure annihilation channel limits for the Galactic Center H.E.S.S. observations are taken 
from [16] and assume an Einasto Milky Way density profile with Qq = 0.389GeVcm^^. Closed contours 
and the marker with error bars show the best-fit cross section and mass from several interpretations of the 
Galactic center excess [17, 18, 19, 20]. 

Figure 2 shows the comparison of the limits from this work with other published limits on 
the DM annihilation cross section. The Pass 8 combined dSph limits are currently among the 
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strongest experimental constraints on the WIMP annihilation cross section and exclude WIMPs 
with the thermal relic cross section for ^ lOOGeV. Our results also constrain DM particles 
with moM above lOOGeV surpassing the best limits from Imaging Atmospheric Cherenkov Tele¬ 
scopes for masses up to ~ ITeV for quark channels and ~300GeV for the T-lepton channel. 

These limits also begin to constrain the preferred region of parameter space for DM interpre¬ 
tations of the Galactic Center excess [19, 17, 20, 18]. However, uncertainties in the Galactic DM 
distribution can significantly enlarge the best-fit regions of (cJv), channel, and widm [21]. 

The future sensitivity to DM annihilation in dSphs will benefit from additional LAT data taking 
and the discovery of new dSphs with current and upcoming optical surveys such as the Dark Energy 
Survey [22] and the Large Synoptic Survey Telescope [23]. The recent discovery of new dSph 
candidates in the first year of the DBS survey [24, 25] suggests that the sample of dSphs that can be 
targeted for gamma-ray searches will continue to grow over the next 5-10 years. Future analyses 
incorporating these new dSph galaxies will improve the sensitivity of the LAT to a DM signal and 
allow the DM annihilation parameter space to be probed to even lower cross sections. 
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